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Abstract 

Different symmetry formalisms for difference eqnations on lattices are 
reviewed and applied to perform symmetry reduction for both linear and 
nonlinear partial difference equations. Both Lie point symmetries and 
generalized symmetries are considered and applied to the discrete heat 
equation and to the integrable discrete time Toda lattice. 

Resume 

Deux formalismes differents pour etudier les symetries des equations 
aux differences finies sur un reseau sont decrits et utilises pour faire la 
reduction par symetrie des equations aux differences finies. Les symetries 
ponctuelles et generalisees sont considerees et appliquees a I’equation de 
la chaleur lineaire discrete et a un treillis de Toda integrable en temps 
discret. 


1 Introduction 


The purpose of this article is to compare two different approaches to the study of 
symmetries of difference equations |p]-p7|,p9 -31 . Both approaches are algebraic, 
in that they both use an infinitesimal formalism in which certain vector fields 
realize a Lie algebra: the symmetry algebra. In one approach the vector fields 
act on the dependent and independent variables and in general on the difference 
equations and the lattices. In the other approach, use is made of evolutionary 
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vector fields, acting only on the dependent variables. In both approaches, further 
choices must be made. In particular, the coefficients of the vector fields can 
depend on the dependent and independent variables at one point of the lattice, 
on a finite number of points, or on an infinite one. They may also depend on 
derivatives of the dependent functions, up to some chosen order, or on finite 
dfferences. 

For point symmetries of differential equations the two approaches are equiv¬ 
alent and the formalisms of ordinary vector fields and evolutionary ones are 
related in a simple manner. In particular, the evolutionary formalism provides 
flows commuting with the original equation. If these flows involve first deriva¬ 
tives only, and they figure linearly, they correspond to point symmetries. The 
two formalisms provide the same symmetry variables and hence the same sym¬ 
metry reductions. 

For difference equations the situation is somewhat different. Also in this 
case, for point symmetries it is very natural to consider vector fields involv¬ 
ing differentiation with respect to both dependent and independent variables. 
Moreover, it is natural to let the transformations act on the difference equation, 
and on the lattice itself. 

A purely difference equation does not contain derivatives (by definition). 
Hence the usual evolutionary formalism, with the characteristic Q of the sym¬ 
metry depending on x, u and derivatives of u, is not a natural tool to use. The 
natural evolutionary formalism for discrete equations is one in which Q depends 
on the values of the variables x and u at different points of the lattice. Such 
evolutionary symmetries are in general not point symmetries and are not related 
to point symmetries in any simple manner. 

One of the most important applications of symmetries of differential equa¬ 
tions is to perform symmetry reduction. For partial differential equations this 
means a reduction of the number of independent variables. 

In this article we shall perform symmetry reduction for various partial dif¬ 
ference equations and show how different types of symmetries lead to different 
results. 

The general theory is discussed in Section 2, first for differential equations, 
then for difference ones. Also in Section 2 we introduce the concept of ’’discrete 
evolutionary vector fields” and their prolongations. Section 3 is devoted to 
symmetry reduction for a linear difference equation, namely the discrete heat 
equation. Translations, as point symmetries, provide a reduced one variable 
equation that is easily solved. Dilations provide a reduction to a dilation-delay 
equation. Still in Section 3 we use discrete evolutionary vector fields to obtain 
solutions of the discrete heat equation, invariant under translations and under 
dilations. Similarly, in Section 4, we use first point symmetries, then discrete 
evolutionary vector fields to obtain reductions of an integrable nonlinear partial 
difference equation: the discrete time Toda lattice. 

2 Continuous Symmetries of Difference Schemes 
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2.1 Equivalence of two symmetry formalisms for differen¬ 
tial equations 

Let us first recapitulate a well-known result for differential equations, namely 
that Lie point symmetries (and also generalized symmetries) can be realized in 
two different ways. The first is by vector fields acting on both independent and 
dependent variables, the second by evolutionary vector fields, acting on the de¬ 
pendent variable only. This is true for arbitrary systems of differential equations 
(ordinary or partial, of any order) [^ j. For point symmetries the two realiza¬ 
tions are entirely equivalent. For generalized symmetries ordinary vector fields 
can always be replaced by evolutionary ones in a straightforward manner. Let 
us illustrate this equivalence on the example of a first order ordinary differential 
equation 

E = Ux - F{u,x) = 0. (1) 

Let us assume that 

X = ^{x,u)dj: + (l){x,u)du (2) 

generates a point symmetry, i.e satisfies 

prXE\E=o = 0. (3) 

The first prolongation of X is 

prX = X + , r=Da:<t>- {D,£,)u, , (4) 

where is a total derivative. 

The prolonged evolutionary vector field, on the other hand, is 

prXe = Qdu + = D^Q, ( 5 ) 

where Q is the characteristic of the vector field. 

The total derivative Dx is itself a generalized symmetry of eq. (|^) (and of 
any equation). Indeed, if we have E = 0, then D^E = 0 follows and so does 
{^Da;)E = 0. Hence, if X satisfies eq. (||) we also have 

{prX - ^D,)E\e=o = 0. (6) 

In view of eq. (^) we have 

prX - + (j)du + [Dx(j) - {Dx^)uj;]du^ - (7) 

“t” E Uxx^Ux] ~ 

= [(t>-^Ux]du + [Dx{(l)-^Ux)]dux- 


With 


Q = (j)-^Ux, ( 8 ) 

we obtain eq (|^) and we have proved that for first order ODEs condition (i is 
equivalent to the condition 

XeE\E=0 = 0 , ( 9 ) 
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i.e. that the two formalisms are equivalent. The same is true for higher order 
ODEs and PDEs of any order [ p^ . 

Given a vector field (^) we can calculate the corresponding characteristic Q 
and the evolutionary field of eq. (^). Hence flows that commute with the 
flow given by the considered equation have the form 

Ux — Q(^Xj Zi, Vjxi '^xx^ 

where A is a group parameter. 

2.2 Point symmetries of difference eqnations 

In two recent articles 0,0 we presented a method for determining Lie point 
symmetries of difference systems. In this terminology a ’’difference system” is 
a system of relations between a set of points in a (p + q) dimensional space 
where the p coordinates xi,...,Xp represent independent variables and the q 
coordinates ui,...,Uq represent the dependent ones. Let us restrict here to the 
case of one difference equation for one function u of one variable x. Let us 
consider a simple case when when three points P,P^,P~ are involved. 

These three points have coordinates (x,u), (x~,u~), and (cc■*■,«+) , respec¬ 
tively. The difference system consists of two relations 

Ea(x,x~,x~^,u,u~,u~*~) = 0, a =1,2. (11) 

This system describes both a difference equation and a lattice. If a continu¬ 
ous limit x~^ X , x~ X exists, then one of these equations goes into a first 
order differential equation, the other into an identity (like 0 = 0). 

The algorithm for finding the Lie point symmetries of such a system is quite 
simple 1^ . We write a vector field as in the continuous case and prolong it 
to all points involved in the difference scheme. In the case of system 0 ) we 
have 

prX = ^{x,u)dx + 4>ix,u)du + ^{x~'~ ,u^)dx++ (12) 

-h (j){x+,u+)d.,,+ +^{x~,u~)dx- + (l){x~,u~)d^-. 

The algorithm for determining the functions ^{x,u) and (j)(x,u) is 

prXEa\Ei=o = 0, 0 = 1,2, i = l,2. (13) 

The algorithm provides functional equations for ^ and </). Solution methods were 
discussed elsewhere, as were applications [0 , |l^ . 

2.3 Commuting flows and evolutionary symmetries for 
difference equations 

Let us again consider a three point ordinary difference system like that of eq. 
(0), but in a form explicitly solved for and x'^: 

El = u'^ — fi{x,x~,u,u~) = 0, E 2 = x^ — f 2 {x, x~,u,u~) = 0. (14) 
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A Lie point symmetry © can be converted into an evolutionary one in the 
same way as was done in the continuous case. Indeed we can define an operator 
Xe and its prolongation prX^ by putting 


prJte = prX - D^+ - ^ D^- 

(15) 

where D^, D^+, and D^- are total derivatives with respect to x, x“*" and x~ , 
respectively. This is equivalent to defining an evolutionary symmetry for differ¬ 
ence equations in the same way as for differential ones, namely 

Xe = [4>{x, U) - ^{X, u)Ux]du 

(16) 

and its prolongation as 


prXe = (</>- Cux)du + - ^~'~u+_^)du+ + 

(17) 


+(</- u^.)du-, 


where the superscripts + and — correspond to total shifts for the corresponding 
variables and functions. Using the operator ( [l^ is equivalent to using opera¬ 
tor (^) for point symmetries. The operator (|l6|) is not a convenient one for 
generalizations going beyond point symmetries. 

Indeed, for difference equations it is more natural (and more fruitful) to 
consider evolutionary vector fields with characteristics Q that depend on the 
independent and dependent variables at different points on the lattice, rather 
than on derivatives, as in eq. & and ©)• In this case we have 


Q{T^x,T^u)du (18) 

prX, = Q{T^x, T^u)du + {TQ)d^+ + {T-^Q)d^- 
m < k < n, m,n € Z. 


Thus Q depends on x and u at a finite, or possibly infinite number of different 
points. In ( |T8|) T is a total shift operator: 

Tx = a;"*", T~^x = x~, Tu{x) = u{x~^) = u'^, T~^u{x) = u{x~) = u~ (19) 

The symmetry algorithm is 

prXeEa\Et=o = 0, a, 6=1,2. (20) 

(together with Eb = 0 we must also use all shifted equations like T^Eb = 0). 
This is equivalent to requesting that the flow 

ux = Q{T'^x,T^u), XX = 0 (21) 

should commute with the flow (©)■ 

To sum up, for difference equations we consider two different algebraic sym¬ 
metry formalisms. The first uses ordinary vector fields of the form (I , coinciding 
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with those used for differential equations. Their prolongation is different, see 
eq. (p^. The algorithm for determining the coefficients of the vector fields is 
given in eq. (liD- These vector ffelds can be integrated to provide genuine point 
transformations taking solutions of difference equations into solutions. The set 
of Lie point symmetries of a difference system is in general much more restricted 
than that of a differential one. This is specially true if we consider a difference 
equation on a fixed lattice, i.e when the lattice equation is just E 2 = — x = a 

where ct is a fixed and nontransforming constant. 

The second formalism for symmetries of difference equations uses the ’’dis¬ 
crete” evolutionary vector fields ( |l8|) They correspond to generalized symmetries 
for difference systems. They provide quite general commuting flows. They can 
be used to perform symmetry reduction. As in the case of differential equations, 
generalized symmetries are particularly useful for identifying systems that are 
integrable, that is those for which a discrete Lax pair exists. For linear differ¬ 
ence equations these ’’discrete” evolutionary vector fields provide commuting 
difference operators 

We demonstrate below, in Sections 3 and 4, that both types of symmetries 
of difference systems are useful for solving difference equations and that they 
provide different types of results. 


3 The Discrete Heat Equation and its Reduc¬ 
tions 

3.1 Evolutionary and point symmetries for the heat eqna- 
tion 

Like the continuous heat equation, the discrete one serves as an excellent exam¬ 
ple of the application of group theoretical techniques. Here we shall use it to 
demonstrate the difference between ’’discrete” evolutionary and point symme¬ 
tries for linear difference equations. 

The discrete heat equation can be written as 


Atu — AxxU = 0 


( 22 ) 


The discrete derivatives used in Ref. were defined as 


Tt-l 




(Xt 


Axx. - 


T" -2Tx + l 


(23) 


( a better notation would have been and A++). Here Tj and are shift 
operators and at and a^ are the steps in the t and x directions, respectively. The 
lattice is fixed, uniform and orthogonal. The Lie algebra of’’discrete” evolution¬ 
ary symmetries was found to be six dimensional (after an infinite dimensional 
subalgebra corresponding to the linear superposition principle was factor out). 
It is isomorphic to that of the continuous heat equation, and was realized by 
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the ’’discrete” evolutionary vector fields 


Po = {Atu)du, Pi = {A^u)du, W = udu, 

B = {2tTf^Aj;U + xT~^u + ]^(JxT~^u)du, 

D = [2tT^^AtU + xT-^A^u + (1 - ]^T-^)u]du, (24) 

K = [t^T^'^AtU + txT^^T~^ A^u + -x'^T~'^u + 


We see that the translations Pq and Pi involve only discrete derivatives, and the 
Galilei transformation _B, dilation D and expansions K involve explicit shifts to 
other points of the lattice. 

Floreanini et. al. obtained an equivalent result for the heat equation 
where At = A^ and A^^ = A~~, i.e. they used ’’down derivatives”, with 


l-T- 


for instance = 

In both cases shifts to a finite number of points are involved. Quite recently 
it was shown [O, |14| that if symmetric discrete derivatives are used, e.g. 


A 


S 

X 


T — T~^ 

-^x 

2(Ja; 


(25) 


then the number of points involved in the symmetries, other than translations, 
will be infinite. Here we will stick with right derivatives, as in ( p^ ) and (24). 

For Lie point symmetries as used in Ref. the definition of discrete 

derivative is immaterial. One simply considers relations between points in the 
lattice. To facilitate the comparison between ’’discrete” evolutionary and Lie 
point symmetries, we shall here consider the following heat equation and lattice 


^m,n+l ^m,n _ ^m+2,n ‘^'^rn+l^n 4“ ^m,n 


,n ,n)^ 

^m+2,n 4“ ^m,n — Oj Xm^n — 0 

^7n,n — O 5 ^m,n+l im,n — ^{Xjn-\-l,n ^m,n) : 


(26) 

(27) 


where c is a constant. The symmetry algebra coincides with the one obtained 
earlier ||l^ using a symmetric derivative on the right hand side of eq. (p^ ) . It is 
spanned by 


Po = dt, Pi=dxj b = xdx + 2tdt, W = udu, S = S{x,t)du (28) 


where S{x,t) is a solution of the system (|26[),(27) and S represents the linear 
superposition principle. The lattice equations (27) can easily be solved and we 
have 


X = (Tjjm + Xq, t = ca^n + to 

where ax, xq and to are integration constants that are not apriori fixed. 
Let us now consider some nontrivial examples. 


(29) 
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3.2 Reductions by Lie point symmetries 


1. Translationally invariant solutions 


A solution of the system invariant under a translation generated 

by Pq — aPi will have the form 


u{x, t) = u{z), z = X + at 

The lattice equations (p^) and the heat equation reduce to 


(30) 


^m+l.n “f ^m — l,n — 0, l,n) (^t) 

— c[zi(Zj77,-|_2,n) 2'u(Zm+l,n) “t” (1^^) 

The solution of eq.(pl|) is 

Zm,n = A(to + Aacn) + zo, (33) 

so Zm,n really depends on just one label N = m+Aacn ( A and zg are integration 
constants). If the reduced equation (|3^ ) is supposed to be a difference equation 
on some lattice, the label N must vary over integer values (N simply enumerates 
different points zg, z±i, z± 2 , independently of their spacing). This implies 
that the constants A, c and a are constrained by the requirement 


Aac = k, k G Z. 

Eq.(|^) can be written as 

u{z + acA^) — u{z) = c[u{z + 2A) — 2u{z + A) + u{z)]. 
The general solution of eq.(p^ can be written in the form 

u(x, t) = cie“^ + C 2 , 

where the constant a is determined by the condition 

- 1 = c[e2“^ - 2e“^ + 1]. 


(34) 


(35) 


(36) 


(37) 


Eq.(p6D also represents a translationally invariant solution of the continuous 
heat equation for a = a. In the continuous limit we have A = (Ta, —> 0 and to 
order eq.(p^) reduces to a = a. 

In general eq.(|3^) is a trascendental equation for a and eq.(p^ determines 
the solution u(z) at a point z + acA^ on the z line in terms of it(z) at three given 
points. As stated above, this will be a point on the same lattice if we choose 
acA^ = kA with k integer. If (Q) is satisfied, then eq.(|^) is an algebraic one 
for V = e°‘^: 


— I = c[v‘^ — 2w + 1]. 


(38) 
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In particular, for fc = 1 eq.(|35D is a three point difference equation and we have 

c+ 1 


aA = In- 


(39) 


2. Reduction by dilation D. 


A scaling invariant solution will have the form 

, _i 

U = U(Z), Z = xt 2 . 

Eq.(|3) imply 


l,n)^ 


Solving eq.(|4l|), we have 


^m,n — 


m — mo 


(40) 


(41) 


(42) 


where toq and uq are constants. E.(^) for Zm,n = const, determines a parabola 
in the (m,n)-plane, so Um,n = u{zm,n) is constant along each parabola. The 
discrete heat equation is reduced to the equation (^^. Choosing some reference 
point z = Zm,n we obtain an equation that can be written as 


/ \ / \ 

u[z -) — u[z) 

In 


c[m(z + 27„) - 2m(z + 7„)+ zi(z)], 7n = —i= , (43) 

^Jc{n - no) 


Even though we are not able to solve eq.(|4^) analytically, we see that a reduction 
has occurred. Eq.(^^ involves one independent variable z, rather than two. 

Eor instance we can take c = 1, rnp = no = 0 in eq.(^^ which then reduces 


to 





(44) 


Eor any fixed value of m we need to give the values of u(z) in a set of equally 
spaced points z, z±7„, z±27„, ... . Eq.(|4^then determines u(z) at irrationally 
spaced points ( ^ along the same line. 

Thus we have reduced to an equation with one independent variable only, 
but it is not a difference equation, rather a difference - delay one. 


3.3 Reduction by discrete evolutionary symmetries 

We rewrite eq.(^) as 


^m,n+l 'am.n — c{VjYn+ 2 .n 4 “ ^771,71)5 C — 




^777+1,77 ^m.n ^3:; ^777+1,77 Im.Tl b; 

^777,77+1 Xm.n — 0, t777^77-|-l ^777,77 — t 


(45) 

(46) 
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and consider reductions of this system by some of the ’’discrete” evolutionary 
symmetries. 

1. Translationally invariant solutions 


The commuting flow corresponding to a general translation is given by (see 
Po and Pi in eq.(|^)): 


= - “ - ^ 

aA (Jf Cx 


(47) 


and an invariant solution will satisfy 


_ / \ 

Wm,n+1 Um,n — ^ V^m+l,n Um^nJ' 
O'x 


(48) 


Together with eq.(^ this implies 


Qj(Um-\-l^n '^m,n) — ('^m+2,n ‘^'^7n-\-l,n T 

This is a linear three point difference equation in m. Its general solution is 

Um,n = A{n) + B{n){l + aax)"^. (50) 

Substituting back into the heat equation (^ ) we find 

A{n + \) = A{n), i?(n + 1) = (1 + a^crt)B(n). (51) 

The general solution of the system (^, hence is 

Um,n = Ci(l + O^UtYll + OCTj,)™ + C2. (52) 

Using the lattice conditions (Q) we can rewrite this as 

u{x, t) = ci(l + acTx)^ (1 + a^at)~ + C2. (53) 


This is not the same solution ( |3^ obtained using translations as point symme¬ 
tries. The two translationally invariant solutions only coincide in the continuous 
limit ax ^ 0, crt ^ 0. 


2. Reduction by dilations 


We shall use the operator D — (1 — ^)W (see eq.(^) to perform the 

reduction. Thus, we solve eq.(|45|) together with the self-similarity condition 


2U 


^m,n ^m,n —1 

i 

o-t 


“h Xr] 


^m.n — 


^ = 0 . 


O'er. 


(54) 
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We solve eq.(|^) for Um,n-i and shift eq.p^) down in n, i.e. replace n by n — 1 
everywhere. Substituting for Um+ 2 ,n-i, Um+i,n-i and Um,n-i we obtain the 
reduced equation 


2c( 71 + l)[Wm+2,n ,n “1“ '^771, n] + m[c(Um+2,n - Um+l,n) - (55) 

2c(Um+l,n “t“ (c + '^ 771 —l,n)] — 0- 

In eq.(|^) m is a variable, n is a parameter. The continuous limit is obtained 
by putting x = Oxfu, t = atn, and taking and at to zero with c = ^ finite. 
We obtain 

2tuxx + xux = 0. (56) 

This is indeed the condition for invariance under dilations generated by D of 
eq. (1^. Eq.([^)is easily solved, eq.(^ is more difficult to deal with. 

We start by integrating eq.(^^ once. The equation involves four values of 
m. We reduce the number to three by putting 


^m,n 


^m+l,n ^m,n 




(57) 


and obtain 

2 c{jl -\- l)(u77j^l 77 T ^[cU 777 _|_l 77 2CU777 77 + (c + 1 )‘^777— 1 ^77] = 0 (58) 

This is a linear ordinary three point difference equation with variable coefficients. 

In order to solve it we use a discrete Fourier transform (often also called the 
Z-tranform), i.e. introduce a generating function 

CX> 

Gn{z) = ^ Vm,nZ"" (59) 

m— — oo 


with 


2TTi 


dzGn{z) 


Cl 2 


m+1 ’ 


where the contour Ci is the unit circle in the complex z plane. Eq. 
implies that Gn{z) satisfies 


(60) 

then 


2 c(n+l)(- l)Gn+[c{Gn z - ) — ‘^CzGn,z + {c+l){z'^Gn,z+zGn)] — 0 . (61) 

z z 


Eq.(^ is easily integrated and we obtain 

(Z-Zi)"(z-Z 2 )’^ 


Gn{z) = 777 


y2n-\-l 


where jn is an integration constant and the complex numbers zi 2 are 

cTz^c 


21.2 = 


c+ 1 


(62) 


(63) 
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Since we have c > 0, z \^2 lie inside the unit circle and eq.(|60|) implies 


— rt ■ T'n 
ZTTl 


dz{z - ziY{z - Z2y' 


Cl 


^ m -\- 2 n -\-2 


(64) 


The dependence on m in eg. (|64|) is explicit, the dependence of 7 „ on n must 
still be determined. To do this we introduce the notation 


Vm.,n = lnlN,n, N = TU + 2n + 2 (65) 

and sustitute into eq.(^^ (the condition for dilational invariance). We also use 
Xm,n = Cxfn, tm,n = CtTi and obtain 

2Tl(Ujji ri Uiii,n—l) T Xfl{Um,n '^m —l,n) — 0- (66) 

In order to introduce Vm,n into eq.(^) we first take the variation (discrete 
derivative) with respect to m and then obtain 

2Tl(Vm,n ^m.n —l) T {jXl T XTlV^fn—l^n — 0* (H^) 

Substituting the expression (^^ for Vm,n we obtain an equation for 7 „, namely 

InliN — l)lN,n — {N — 2n— 2)/7V-l,n] ” N-2,n-l = 0. (68) 

Eq.(^^must hold for all values of N and 7 „ is independent of N. The quantity 
lN,n is defined by the integral in eq.(|^) and can be evaluated using the residue 
theorem. For general N and m this is not very illuminating. For low values of 
N and n > 0 the results are quite simple, for instance 


^l.n — 

{ziz2r = i^r, 

^Q.n — 

0, 


0, 

l2,n = 


II 

e 


II 

e 

-|n(n-l)(^)"-i(Hll: 


Substituting these values into eq.(68) we obtain 


7„ = (c+ l)"7o, 


(69) 


which can be shown to be valid for all integer values of n, both positive and 
negative. Finally, the m variation Vm,n of the dilationally invariant solution of 
the discrete heat equation (see eq.(^) satisfies 


Xm,n — bo(o T 1) ^ . 

ZTTl 


dz{z - zi)^{z - Z2y 

N 


/Cl ^ 

where Ci is the unit circle and zi ,2 are defined in eq.( 


(70) 
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The continuous limit of this self similar solution is 


v{x,t) = —u(x,t) = (71) 

dx 

which satisfies the continuous limit of eq.(|54[), namely 2tut + xux = 0, which 
implies 

2tVx + XV = 0. (72) 

4 Symmetry Reduction for the Discrete Toda 
lattice 

A recent article was devoted to a hierarchy of nonlinear integrable difference 
equations associated with a discrete Schrodinger spectral problem [^. The 
hierarchy includes the discrete Toda and discrete Volterra lattice equations. All 
equations in the hierarchy involve two independent variables: discrete space and 
discrete time. 

In this section we shall use Lie point symmetries and ” discrete” evolutionary 
symmetries to perform symmetry reduction for the simplest equation in the 
hierarchy, namely the Discrete Time Toda Lattice (DTTL) itself: 

g«n,m-Mn,m + l _ + l-“n,m + 2 _ q, 2 ^gUn-1 ,m + 2-“n,m + l _ g«n ,m+l ““n + l ,m ^ 

( Of is a constant). 


4.1 Reduction by Lie point symmetries 

We complement the DTTL (|7^ ) by the lattice equations: 


^n+l,m H” ^n —l,m 

= 0, 

^n,7n — ^n,m+l5 

(74) 

^n,m+l H” ^n,m—1 

= 0, 

^n,m — 


The solution of eq.(Q is 




Xn,m — 

in,m 

= (Jtm, 

(75) 


where the integration constants and at represent the lattice spacings and we 
have set the two further ’’initial value” constants xo.m and tn,o equal to zero. 

The Lie algebra of the Lie point symmetry group of the difference system 
( [ 73 I ) and ( [ 7 ^ ) is spanned by 

Po = dt, Pi = dx, Do = tdt, Di = xdx, W = du- (76) 

Let us look separately at reductions by translations and by dilations. 

A. Translationally invariant solutions. 
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A general translationally invariant solution has the form 

'^n,m = ^n,m) — 

^n,7n — ^n,m “t“ Ciin^m — ^x'^ “1“ (XCTiTTl. 

Eq.(|^) reduces to an equation for it(^) evaluated at five points namely 

^n,m — ^7 Cn,m+1 = ^ “t“ Q’O't^ ^n,m+2 — ^ “1” (^8) 

^n+l,m — C A ^X 7 Cn — l,m+2 — C “t“ ‘2ciO't- 

As in the case of the heat equation (see eq. ( p4p^ )) the reduced equation will 
be an ordinary difference equation on a lattice only if ^n,m is a function of one 
discrete label. This occurs if we put 

aat = k(Tx, k G Z (79) 

in eq. In particular, if we choose A: = 1 in eq.(|7^) we obtain a three point 
difference equation (since ^n.m+i = ^n+i.m = Cn-i,m+ 2 )- The right hand side 
of eq. dzl) vanishes identically and the left hand side implies 

‘^n,m ‘^'^n,m+l “t“ 1ln,m+2 — 0- (^0) 

The general solution of eq.(|80|) is 

Un,m = f{n)m + g{n). (81) 

Substituting back into eq.(|7^) we obtain the general translationally invariant 
solution of eq.(|7^) as 

Un,m = An{n + m) + Bm + Cn + D, (82) 

where A, B, C and D are constants. 


B. Reduction by dilations 

A solution invariant under dilations generated by Dq — (iDi will have the 
form 

’^{^71^1717 ^n,m') — Cn,m — = C' (^^) 

Substituting into the DTTL ( [^ we obtain a nonlinear equation involving u(^) 
evaluated at 5 points: 

with ^ as in (|^) and 

{.={(^)t & = (85) 

m m 

6 =C( ”^'^^ )^-o'a:Crf(TO + 2)^,^4 = ^ + 

m 

Eq. (Q) is a dilation - delay equation. 
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4.2 Discrete Evolutionary Symmetries for the DTTL 

The DTTL can be written as a system of two equations as 

( 86 ) 

(87) 

where we have 


^n,m+l ^n,m — Q^(^n,m+1 ^n+l,m) 


^n,m+l 


^n,m — 


'^n^m — ^n,m’^n+l,m5 '^n,m — ^ Ti,,m^ 

TT — n • 


( 88 ) 


As in the case of the Toda lattice itself, both isospectral and nonisospectral 
’’discrete” evolutionary symmetries exist ^ and we shall consider both 

of them. 

The simplest nontrivial isospectral (generalized) symmetry is given by 


(Un,m)e — Un,??! [Un—l,m Un+l,m ~t“ ^ri,m ^n+l,m] (^^) 

(bn,m)e — 0 -n—l,m[bn,m + bn—l,m] Cln,m[bn +1 ,771 H“ ^n,m] ■ 

We assume a„,m 7 ^ 0. Symmetry reduction amounts to solving eq.(^^ for 
{an,m)e = {bn,m)e = 0. Eq. (^ Can be once integrated to yield 


^n—l,m T ^n,m “b ^n,m — 
an,m(^n+l,m “b ^n,m) — 


(90) 


We can eliminate bn^m from (|^) to obtain 


I ( \l^n,m ^n+l,7n “b \l-^m ^n — l,m ^Ti,m} — 


(91) 


Eq. (0) can be viewed as a discrete analog of the equation for elliptic functions. 
It is an integrable difference equation in one variable n with m as a fixed param- 
ether. The functions Am and Bm must be determined by putting the solution 
ofeq.(|§) into (@. 

The simplest nonisospectral symmetry is given by 


(an,m)e — ‘2l7l (277. “b 2777 l)^n,m] (^^) 

(^n,m)e — 4 T 2[(77 Tfl -\- l)77n,m “b 777 l)Uj^_i^m], (Ull) 

( see eq. (64) in Ref. Q with A: = 0, a = 1 ). 

Again, we must solve the equations {an,m)e = {bn,m)e = 0. We eliminate 
bn,m from these two equations and obtain 

(277 -\- ‘2tTI 3)^(77 “b 777 “b “b (114) 

— [(77 + m){2n + 2m + 3)^ + (77 + 777 + 1)(277 + 2m — l)^]a„,m + 

+(77 + 777 — 1)(277 + 2m — l)^a„_i^m = 16(277 + 2m + 1). 
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One solution of the corresponding homogeneous linear ordinary difference 
equation (in which n is the independent variable, m a parameter) can be guessed, 
namely 

(95) 


1 


{n + m){n + m + 1) 

The general solution of the inhomogeneous equation is constructed in the form 

0‘n,m — ^71,711^1^,711, ■ (96) 

We define 

Cn,m — Pn+l ,m - Pn ,m ■ (97) 

Substituting (Dm into (p® we obtain a first order equation for c„_m, namely 


(n + m + 2)(2n + 2 to + 1)^ ^ ^ + m + 2){2n + 2m + 3) 

Cn+l,m — 7 7T77; ~ ZT^^Ti.m T 16- 


(n + m + l)(2n + 2 to + 5 )^ 


(2n + 2 to + 5)2 


(98) 


A particular solution of the homogeneous part of eq. (|98| ) is 

Cn,m = 2(n + m + 1) (99) 

which implies = n{n + m + 1). Finally, the general solution of eq.(p^ is 
1 r w , B{m) 


^n.m — 


(n + m)(n + m + 1) 


{A(m) 


{2n + 2m + 1)^ 


+ n(n-I-2m-|-1)} (100) 


where A{m) and B{m) are integration constants. Returning to (93) for (6„,rn)e = 
0 we find 

( 101 ) 


(2n + 2 m-l)(2n + 2m+l)' 


invariance 


Together, eq. and ( |101[ ) represent the general solution of the ir 
condition (a„,m)e = (&n.m)e = 0 in eq.(||, 

The functions A(m) and B(m) must be determined by substituting Un^m and 
5n,m into the DTTL of eq. (^) (for 0=1). This can be done; the result is: 


B{m)=0, A(m) = m(m + 1) 


( 102 ) 


and finally 

^ 71,771 — I; ^ 71,771 — 0. (103) 

The result is somewhat disappointing since the invariant solution is one in which 
the fields are equal to their asymptotic values for all n and m. 

In order to obtain more interesting solutions, higher symmetries must be 
considered. Those will however be nonlocal and that is beyond the scope of the 
present article. 
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5 Conclusions 


We have shown how to use symmetries of difference equations in two variables, 
either Lie point or ’’discrete” evolutionary ones, to construct solutions by car¬ 
rying out a symmetry reduction. 

We are confronted with two different situations. In the case of Lie point 
symmetries the symmetry variables, as opposed to the case of partial differential 
equations, do not always reduce the equation to a difference equation in a space 
of lower dimension. We have to impose a further constraint to be able to do so. 
Moreover in some cases, among them that of dilations, the system reduces to a 
dilation - delay equation which is difficult to solve. 

The situation is different for ’’discrete” evolutionary symmetries. These 
symmetries exist only when the discrete system is integrable and the symmetries 
do not act on the lattice. The symmetry reduction can always be carried out, 
but the obtained reduced equation can be nonlinear and difficult to solve (this 
is also true in the continuous case). 
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